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Abstract 

Various experimental settings that involve drying solutions or suspensions of nanoparticles - often called 
nanofluids - have recently been used to produce structured nanoparticle layers. In addition to the for- 
mation of polygonal networks and spinodal-like patterns, the occurrence of branched structures has been 
reported. After reviewing the experimental results we use a modified version of the Monte Carlo model 
first introduced by Rabani et al. [Nature 426, 271 (2003)] to study structure formation in evaporating films 
of nanoparticle solutions for the case that all structuring is driven by the interplay of evaporating solvent 
and diffusing nanoparticles. After introducing the model and its general behavior we focus on receding 
dewetting fronts which are initially straight but develop a transverse fingering instability. We analyze the 
dependence of the characteristics of the resulting branching patterns on the driving chemical potential, the 
mobility and concentration of the nanoparticles, and the interaction strength between liquid and nanoparti- 
cles. This allows us to understand the underlying instability mechanism. 
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I. INTRODUCTION 



Branched structures are ubiquous in nature. They appear in a wide variety of phenomena in- 
cluding growing trees, river networks, dendrites in solidification, and crystallization Q. In hydro- 
dynamical systems they occur, for instance, in Saffmann-Taylor fingering [2] whose zero-interface 
tension limit is well described by the discrete computational model of diffusion limited aggrega- 
tion (DLA) [3 j. Branched patterns may also result from contact line instabilities in spreading drops 
of surfactant solution Hill 16). This is, however, rather an exception as in most cases transverse 
instabilities of advancing or receding contact lines do not lead to branched structures. Normally 
one finds either (i) arrays of straight parallel or wedge shaped advancing fingers, as for a gravita- 
tionally driven front moving down an incline or vertical wall [H |9l QUI [H) or driven by thermal 
gradients ll2l[T3l[l4]|; (ii) fingers which advance radially outwards in the case of a spinning drop 
IfBH : or (iii) fingers (Tgl El [H or fields of droplets flU [J9l EOl ED remaining at rest behind a 
receding circular or straight dewetting front (cf. also section 3 of Ref. Il22l0 . 

For the presently studied system of an evaporating solution of gold nanoparticles that dewets 
a silicon substrate (23l|24l[25l[26l[27l|, fingering was predicted Il28ll using a kinetic Monte Carlo 
approach [[26l [29l [30). Recent fine-tuned experiments reviewed below in Section [II] have found 
such instabilities, analyzed their dependence on the properties of the nanoparticles, and could 
furthermore justify the usage of the seemingly simplistic Monte Carlo model II3TT1 . This will be 
elucidated in more detail below in Section HITAl 

In the broader field, dewetting experiments with evaporating suspensions of particles or 
macromolecules are performed by various groups in a wide variety of settings. Examples in- 
clude drying macroscopic drops of colloidal solutions of micron- or nanometer- sized particles 
J32l[33l[34l|35l|36l, spin-casted aqueous collagen solutions J37l|38l|39]], spin-casted solutions of 
polystyrene in benzene [40J, and spin-casted solutions of poly acrylic acid (PAA) in toluene PTTl . 
The behavior presented in these studies is generic for a wide class of solvents and solutes (see, 
e.g., references in [|40l0 . More complex situations have also been studied, such as evaporating thin 
films of a binary polymer solution on horizontal ll42ll or inclined [43] substrates, a drying binary 
suspension of hard- sphere colloidal particles and a non- adsorbing polymer fi4ll . and structuring 
of polymer solution films caused by high-temperature evaporation/boiling in a dip-coating setting 

For the less complex systems, in several regions of the parameter space, the dewetting behavior 
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of the solutions appears to be very similar to that described for simple polymeric liquids [[22l |47l 
58l|49l. In these regions it is possible to use the solute to 'image' the dewetting patterns created 
by the volatile solvent. This route was taken, e.g., in Refs. 071 [39]]. This is, however, not the case 
in other parameter regions. Most strikingly, the experiments show that dewetting (evaporating) 
solutions show much stronger transverse instabilities of the dewetting front than pure liquids (see, 
e.g., Refs. EH gD El) 

The present paper focuses, after a short review of the experiments (Section [II]), on an analysis 
of the fingering instability employing the kinetic Monte Carlo model as developed in Ref . [29 ] and 



later modified in E6ll . Section [HI] introduces the model and the numerical algorithm. Section IV 



reviews the general behavior of the model, whereas Section [V] entirely focuses on a parametric 



analysis of the fingering instability. Finally, Section [VI] concludes and gives an outlook of future 
work. 



II. EXPERIMENTS 

In the following, we review recent experiments that use monodisperse colloidal solutions of 
thiol-passivated gold nanoparticles, sometimes called a 'nanofluid'. The 2 or 3nm gold core is 
covered by alkyl-thiol molecules, whose carbon chain length can be changed from C 6 to Ci 2 EE). 
Varying the length of the chain allows control of the interaction parameters, to a certain extent. 
In general, the particles are hydrophobic and are then easily suspended in toluene. They can also 
be made hydrophilic by modifying the termination of the passivating alkyl-thiol molecules. Then 
they can be suspended in water or in another polar solvent. Most experiments are carried out using 
hydrophobic particles. Normally, the solution is deposited (see below) on a silicon substrate that 
is only covered by the native silicon oxide layer E51 . However the wetting behavior of the solvent 
(and the interaction of the nanoparticles with the substrate) can be (locally) changed by oxidizing 
the substrate further ll27ll . The properties of the solvent can be changed as well, e.g., by adding 
excess thiol (which is thought to mainly affect the viscosity of the solution lt3TT0 . 

The actual deposition process is very important to determine evaporation speed and, in con- 
sequence, the results of the ongoing pattern formation process. Two procedures are followed. 
On the one hand, a thin film of solution is deposited by spin-coating a drop of solution onto the 
native-oxide terminated silicon substrate. Then evaporation competes with dewetting as was de- 
scribed some time ago in Ref. [39] for a related system involving a solution of macromolecules. 
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After all solvent has evaporated the resulting nanoparticle deposits are imaged using atomic force 
microscopy (AFM). Depending on the concentration of the solution one may find mono-modal 
or bi-modal cellular network structures as in Ref. [25], ribbon-like or labyrinthine structures, or 
branched structures (see Fig.[T]). Flower-like branched structures are also observed (Fig. [2]). Note 
that for spin-coated films the evaporation is fast and normally the structuring is complete even 
before the spin-coater is stopped. At present, structuring of the film under spin-coating conditions 
has not been observed in situ. 




FIG. 1: AFM images of dodecanthiol-passivated 2nm gold nanocrystals, spin-coated from toluene onto 
native oxide terminated silicon substrates. From top left to bottom right, concentrations increase as follows: 
O.lOmg/ml, 0.25mg/ml, 0.50mg/ml, 0.75mg/ml, l.OOmg/ml, 1.25mg/ml, 1.50mg/ml and 1.75mg/ml. Scans 
are 2 x 2 to 5 x 5^m 2 in size. All scale bars correspond to lfim. The empty substrate is white, the deposited 
nanoparticles are black. 

A meniscus technique using a teflon ring was recently introduced II5TT1 to decrease the solvent 
evaporation rate during nanoparticle self-organization. Precursors of this technique have involved 
the use of latex beads (see, e.g., Ref. Il52l and references in Ref. ll5TT0 . The slowing-down of the 
evaporation process improves control and allows the use of contrast-enhanced microscopy ll53ll to 
capture the dewetting process on video II3T11 . In the meniscus technique a drop of solution is de- 
posited onto a teflon ring that sits on the silicon substrate. As the toluene wets teflon, it evaporates 
quickly (within seconds) from the center of the ring only. The meniscus thereby formed between 
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FIG. 2: AFM pictures of flower-like branched structures in the case of (a) a toluene solution of gold nanopar- 
ticles spin-coated onto native oxide terminated silicon substrates, and (b) a toluene solution (with an excess 
of thiol) deposited using the meniscus technique. The empty substrate is white, the deposited nanoparticles 
are black. Scale bars correspond to 2/im. 

the silicon substrate and the teflon ring slowly shrinks due to evaporation on the time scale of one 
hour. The pattern formation process is confined to the region of the receding toluene/silicon/air 
contact line. 

Using the meniscus technique one can observe (as for the spin-coating procedure) labyrinthine 
spinodal structures and mono-modal or bi-modal network patterns depending on nanoparticle con- 
centration and position inside the teflon ring. (The position within the ring is closely linked to the 
solvent evaporation rate). It therefore becomes possible to study branched structures in a more 
controlled manner. Selected examples of such structures are given in Fig. [3} The systematic study 
performed in Ref. PTTl indicates two important effects: Fingering strongly depends (i) on the chain 
length of the thiol molecules passivating the gold core of the particles and it also depends (ii) on the 
amount of excess thiol in the solution. The chain length is varied from pentane (C 5 ) to tetradecane 
(C14) at fixed nanoparticle concentration and deposited volume [see Figs. 2(a)-(e) of Ref. 1(3111 1. 
For C 5 and C 8 ligands, no formation of branched structures is observed irrespective of evaporation 
rate, i.e., position within the teflon ring. For longer chains (C10 and C12) well-developed branched 
structures are formed. Increasing the chain length even further (C14), however, produces less de- 
veloped branching. An addition of a small amount of excess thiol to the solution (0.1% by volume) 
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FIG. 3: AFM pictures of fingering patterns obtained by evaporative dewetting using the meniscus technique. 
The employed thiol-passivated gold nanoparticles are coated with thiol molecules of different lengths. 
Shown are images representative of chain length (a) Cg, (b) Cio, (c) C12, and (d) C14. No excess thiol 
is present. The empty substrate is white, the deposited nanoparticles are black. Scale bars correspond to 
1/xm. 

enhances the development of branched structures strongly [see Ref. ED]]. For small chain length' 
(C5- and C 8 ) branching (albeit sometimes weak) is now observed. The branching is considerably 
stronger for Cio, C i2 and C14. 

The main aim of the present work is to shed more light on the mechanisms underlying the 
branching process using a simple computational model J26l [29j. To separate the process from 
other effects we Analise the dependence of the branching of a straight evaporative dewetting front 
on the control parameters. The results are used to explain the experimental findings and indicate 
avenues for future experiments. We will come back to this below in the Conclusion. 
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III. KINETIC MONTE CARLO MODEL 



A. Justification of usage 

A simple two-dimensional kinetic Monte Carlo model for the present process was first pro- 
posed in Ref. Il29l to model experiments such as those presented in Ref. G4) . Before we introduce 
the model in detail, we first discus its applicability and justify its usage. The model is based on 
two key assumptions: (1) that the relevant processes can be captured by a two-dimensional setting 
neglecting changes in the film thickness of the evaporating film; and (2) that all relevant dynamics 
are a result of the diffusion of nanoparticles and of the evaporation of the solvent. Convective 
motion of the solution is entirely neglected. A refined model increases the spatial range of con- 
sidered energetic interactions including next-nearest neighbors E6l but otherwise uses the same 
basic assumptions. Considering the strong assumptions, the agreement with experiments is amaz- 
ingly good J26l [29]]. As we shall now discuss, recent experiments using the meniscus technique 
(described above) may explain why. 

In Ref. OTTl the evolution of the branched patterns is followed in real-time using contrast- 
enhanced video microscopy. The video (complementary material of Ref. [31]) clearly shows that 
different processes occur on different scales. First, a macroscopic dewetting front recedes, leaving 
behind a seemingly dry substrate. The macroscopic front can be transversally unstable resulting 
in large-scale (> 100/zm) strongly anisotropic finger structures. For fronts that move relatively 
quickly these macroscopic structures cover all available substrate. However, when at a later stage 
the macroscopic front becomes slower, those fingers become scarce and 'macroscopic fingering' 
finally ceases. At this stage it is possible to appreciate that the seemingly dry region left behind by 
the front is not at all dry, but covered by an ultrathin 'postcursor' film that is itself not stable. At 
a certain distance from the macroscopic front the ultrathin film starts to evolve a locally isotropic 
pattern of holes. The holes themselves grow isotropically in an unstable manner resulting in an 
array of isotropically branched structures as shown, e.g., above in Figs. [2] and [3] This indicates that 
nearly all of the patterns described in various publications result from processes in the ultrathin 
film whose thickness is of the order of the size of the nanoparticles. 

Focusing on the structuring of the ultrathin film, we next identify the important dynamical 
processes. On the mesoscopic scale the dynamics of a thin film of pure liquid can be described 
by a thin film equation derived using a long wave approximation [|54l [55). In such a model, the 
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temporal change in film thickness h results from the gradient of the convective flow h 3 (Vp)/3r] 
and from an evaporative loss /3(ppo — p)/p [HH [57]]. Here, r] and p are the dynamic viscosity 
and density of the solvent, respectively, /x is the chemical potential (related to the ambient vapor 
pressure of the solvent), and (5 is a rate constant that can be obtained from gas kinetic theory or 
from experiment 1551 . The pressure p contains curvature and disjoining pressures and drives both 
processes - convection and evaporation. The disjoining pressure describes the wettability of the 
substrate [J22l |47l [58l [59J. Due to the thickness dependence /i 3 , the mobility related to convective 
motion is large for thick films but decreases strongly with decreasing film thickness whereas the 
mobility related to evaporation (the rate constant) is a constant. In consequence, it is expected 
that for a nanometric film the evaporation term strongly dominates. Using the parameters given in 
Ref. Il56ll and assuming a small contact angle of about 0.01, one obtains a cross-over thickness in 
the lower single digit nanometer range. Below this thickness the solvent dynamics is dominated 
by evaporation. 

This consideration, together with the experimental observation discussed above, justifies the 
neglect of convective motion in the kinetic Monte Carlo model. Without convection a two- 
dimensional model is sufficient to cover the essential processes. Note that a refinement proposed in 
ll27l introduces a three-dimensional aspect into the two-dimensional model by making the chem- 
ical potential dependent on mean liquid coverage (i.e., on a parameter related to mean thickness). 
This amounts to a consideration of the thickness-dependent disjoining pressure in the evaporation 
term without the explicit incorporation of a film thickness. The resulting pseudo three-dimensional 
model successfully reproduces bi-modal structures (23. These shall, however, be of no concern 
here. 

B. The model 

Kinetic Monte Carlo models are a common tool to investigate a wide variety of dynamical 
processes j60l (6TJ. In particular, they have proven to be attractive tools to model adsorption, 
diffusion, and aggregation in surface science resulting in interfacial growth and structuring like, 
e.g., the growth of metal layers as atoms arrange after being deposited on surfaces Il62ll63ll64ll . 

The approach followed in Ref. E9l to model an evaporating dewetting nanoparticle solution 
is based on an Ising-type model for the liquid-gas phase transition. To facilitate comparison with 
previous work we here choose the model to be exactly the same as in Ref. E6ll (building on 
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FIG. 4: Schematic picture of the two-dimensional lattice gas model used: the size of nanoparticle is in our 
case 3x3, whereas a liquid cell is of size lxl. The side length of such a cell represents the correlation 
length of the solvent approximated to be 1 nm ll29l . 

Ref. Il29l0 : The system is described using a two-dimensional lattice gas of two fields n and /. We 
allow for three types of possible states of a cell - liquid (Z = 1, n = 0), nanoparticle (/ = 0, n = 1), 
and vapor (/ = 0, n = 0, i.e., cell empty). Thereby liquid/vapor is assigned per lxl cell whereas 
it is assumed that the nanoparticles fill 3 x 3 cells (see sketch Fig. [4]) to indicate the larger size of 
the nanoparticles used in the experiment as compared to the correlation length of the liquid. This 
detail, however, turns out to be not very relevant as further discussed in the Conclusion. 

The energy of a configuration of the two-dimensional assembly of vapor, liquid and nanoparti- 
cles is determined by the Hamiltonian 

E = J2 n . n . - S -jY^ n * l J ~ y Yl ~ (!) 

<ij> <ij> <ij> i 

where E\u e nn and e n \ are the interaction energies for adjacent sites filled by (liquid,liquid), 
(nanoparticle,nanoparticle) and (liquid,nanoparticle), respectively. For fixed interaction strengths, 
the equilibrium state is determined by the chemical potential /i that controls the evapora- 
tion/condensation of the liquid. We will call the driving force as it is primarily responsible 
for the motion of a dewetting or wetting front. The sums J2<ij> are taken over all pairs of nearest 
and next-nearest neighbors. Thereby, the interaction strength of the next nearest neighbors is cor- 
rected by a factor l/y/2 due to their larger distance G6ll . In the following we fix su = 1, i.e., we 



express all energies and the chemical potential in the scale of the liquid-liquid interaction energy. 

The energy functional determines the equilibrium state and the energy landscape of the system. 
The dynamics is determined by the allowed Monte Carlo moves, their relative frequency, and the 
rules for their acceptance. Two types of moves are allowed: (i) evaporation/condensation of liquid 
and (ii) diffusion of nanoparticles within the liquid. To simulate the dynamics, the energy loss or 
gain AE related to a potential move is calculated. The move is then accepted with the probability 
Pace = min[l, exp(— AE/kT)] where k is the Boltzmann constant and T the temperature. The 
rule implies that any move that decreases the energy of the system is accepted when tried. Any 
move that increases the energy has a probability corresponding to the related Boltzmann factor. 
Consistent with the above scales, temperature is expressed in units of en = 1/k. It determines 
the importance of fluctuations in the system. For T = the system is fluctuation-free, i.e. the 
evolution follows a deterministic gradient dynamics. 

Practically, we use a checkerboard Metropolis algorithm to advance the liquid/vapor subsystem; 
each solvent or vapor cell is examined in turn, and is converted from liquid to vapor or from vapor 
to liquid with the acceptance probability p acc . After one solvent cycle all particles are considered 
and a diffusive move is tried. The nanoparticles are free to perform a restricted random walk on the 
lattice. It is restricted because particles are only allowed to move into wet areas of the substrate, 
i.e., onto cells with 1 = 1. This is a very important rule and models zero diffusivity of the particles 
on a dry substrate. If such a 'forbidden' move is tried, the particle is left at rest and the next particle 
is considered. The diffusivity or mobility of the nanoparticles inside the liquid is controlled by the 
number of times each particle is examined after one solvent cycle. This computational ratio M of 
particle and solvent cycles reflects the physical ratio of time scales for evaporation and diffusion. 
Large M, for instance, indicates that diffusion is fast as compared to evaporation, i.e., it stands 
for a large diffusion constant of the nanoparticles or/and a low evaporation rate for the liquid. It 
therefore represents a low effective viscosity of the solution. 

The model is used in the following to first review some general results for a homogeneous 



system without and with nanoparticles (Section |TVj). Then we focus on a parametric analysis of 
the fingering instability (Section fV|). 
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IV. BEHAVIOR OF A HOMOGENEOUS SYSTEM 



A. Without nanoparticles 



Without nanoparticles the behavior of the model introduced in Section III is well known, it 
reduces to the classical Ising model in two dimensions describing the order-disorder transition in 
a ferromagnetic system. Using the mapping /i — > [i H — 2 and I — > (s + l)/2 where s = ±1/2 
corresponds to the spin of a cell, /i to the magnetic moment and if to the external magnetic field, 
all results since Onsager ll65ll can be used. In particular, an infinitely extended system has a critical 
point at \l c = -2 and /cT c °° = l/[2 ln(l + y/2)] « 0.567. For temperatures below T c °° at /x ph = -2 
the liquid state and the vapor state may coexist, whereas for /i < — 2 [/i > —2] eventually the 
vapor [liquid] state dominates. A chemical potential not equal to -2 corresponds to a non-zero 
external magnetic field in the ferromagnetic system. 

For T < T c °° there exists a first order phase transition at /x ph ; systems in the liquid [gas] state are 
metastable for a certain range below [above] /x ph . Mean field theory as presented, e.g., in Ref. Il66ll 
provides the line that is the lower [upper] limit for the existence of metastable liquid [gas] states 
in the (/x, kT) plane. The two curves fulfill 

^ = i4TU-l -2, (2) 



ms — g l rp 

where we set the grid spacing used in Ref. [66] to one. 

An overview of the equilibrium phases and the resulting behavior of a straight wet- 
ting/dewetting front (as described next) is given in Fig. [5} Starting with a liquid covered substrate 
in a region (/i < — 2) where the global energy minimum corresponds to the vapor phase leads to 
an evaporative dewetting process that either follows a nucleation and growth pathway or results 
from a spinodal-like process. Typical snapshots illustrating the two processes are shown in Fig. [6} 
Nucleation and growth of holes occurs in the parameter region close to /x ph = — 2 where the ho- 
mogeneous liquid state is still metastable. At smaller /i < /i~ s , i.e., at larger driving forces 
a spinodal-like process occurs when starting with a liquid-filled plane. Many very small holes 
appear at once, and all liquid evaporates very quickly, leaving the substrate empty. When starting 
with an homogeneous gas state at /i > — 2, similar mechanisms result in a liquid-filled plane. The 
relevant border there is \i = /x+ s . 

A straight front separating vapor- and liquid-covered substrate areas remains on average at rest 
for /x = —2 although it might fluctuate quite strongly. The liquid state will advance and recede 
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FIG. 5: Phase diagram without nanoparticles (central panel). A schematic overview of system behavior 
is given in the plane spanned by the dimensionless chemical potential /i and the dimensionless measure 
of temperature kT. The liquid and gas phase are indistinguishable above the critical temperature kT c = 
l/[21n(l + \/2)] « 0.57. The corresponding gradual change of mean density with /i is indicated by the 
four snapshots of equilibrium states shown in the upper line (from left to right: /i = —2.75, —2.25, —1.75 
and —1.25; kT = 0.8). Below T c the system can show a rich dynamics when evolving towards equilibrium. 
For \i > —2 [fi < —2] the equilibrium corresponds to liquid [gas] as illustrated by the snapshots in the 
bottom row (fi = —1.5 [fi = —2.5] and kT = 0.05), whereas at /x p h = —2 gas and liquid may coexist. At 
/Xph a straight front separating gas and liquid does not move on average. However, depending on temperature 
it fluctuates, as illustrated by the snapshots in the right column (from bottom to top kT = 0.1, 0.3, 0.45 and 
0.55; ji = /x p h = —2). A straight liquid front will recede [advance] for ji < — 2 [fi > —2], i.e. one is able 
to study evaporative dewetting [wetting] fronts. Starting, however, with an homogeneous liquid-covered 
[gas-covered] substrate for /i < — 2 [fi > — 2] the substrate will empty [fill] via a nucleation or spinodal- 
like process. The borders between the different processes are indicated as heavy lines and are discussed in 
the main text. All snapshots are obtained from simulations for a small domain size of 300 x 300 after 1000 
kinetic Monte Carlo steps. 
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FIG. 6: Snapshots from typical evaporative dewetting processes without nanoparticles in the case of (a) a 
spinodal-like process at /i = —2.55 and (b) of nucleation and growth of holes at \i = —2.3. The starting 
condition in each case is an homogeneous liquid film. The number of Monte Carlo steps is from left to right 
(a): 5, 10, and 15; (b): 40, 80, and 120; kT = 0.3 and the lattice size is 600 x 600. Liquid is gray (green 
online) and the empty substrate is white. 

for fi > —2 and \i < — 2, respectively. For a resting or moving front the temperature (T < T c ) 
determines the strength of the fluctuations that modulate the straight front (see snapshots at the 
right of Fig. [5]). The overall picture is similar for circular fronts, however, the exact value of /i 
that results in a resting front depends on the average curvature of the front (i.e. the size of the 
hole). A two-dimensional 'curvature pressure' enters the balance. Fixing some /x~ s < /x < —2, a 
hole in a liquid layer will grow [shrink] above [below] a certain radius r c (/x). The growing hole 
will remain almost circular. Practically, on a square grid in late stages it will become square-like 
(original model in Ref. [29] where only nearest neighbors are considered). Including next nearest 
neighbors leaves the holes circular up to a later stage, when they become octagonal. 
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B. With nanopar tides 



Having discussed the behavior of the model without particles we now turn our attention towards 
situations where nanoparticles are present. First, we discuss the equilibrium behavior, followed by 
an analysis of the evolution of an initially liquid-filled well-mixed homogeneous system. 

We expect the particles to influence the location of the liquid-gas phase transition. Inspecting 
the Hamiltonian ([T]) one can estimate the influence using a mean field argument. Let us consider 
a liquid-filled cell at a straight liquid front. We consider the energy needed to move the front, i.e., 
to empty the cell. On average a cell at a straight front has two direct neighbors occupied either 
by liquid or nanoparticles. We replace the nanoparticle occupation number in the liquid-particle 
interaction term by its mean value, i.e., the coverage 0, neglect the particle-particle interaction 
term as we consider a cell either filled by liquid or gas, and replace one of the liquid occupation 
numbers in the liquid-liquid interaction term by its mean value for a filled cell 1 — (f). In this way 
we can map the model onto the pure liquid-gas case by replacing fi by jl = /x + 2{e n i — 1)0. 
Following this argument we expect the phase transition to occur at jl = —2, i.e., an attractive 
liquid-gas interaction (e n i > 0) stronger than the liquid-liquid interaction (su = 1) implies that 
the phase transition occurs at /2 ph = — 2 — (e n i — l)cp < /i ph . Note, however, that for low particle 
concentrations (e.g., = 0.1) the change is rather small for e n \ = 1.5, i.e. the value used in most 
of the present work. This agrees with simulations where we have not spotted a significant shift for 
the used concentrations < 0.2. 

The aspect that most interests us here is, however, not the equilibrium behavior but the dy- 
namics of the liquid-gas phase transition. As in Section [TV A| we start at a /i < —2 with a liquid 
covered substrate, i.e., in a region where the final state without nanoparticles corresponds to the 
vapor phase. Now, however, although the liquid evaporates the nanoparticles remain. Depend- 
ing on the particular parameter values chosen one finds final structures that do not change any 
more on short time scales. They might, however, slowly evolve, e.g. by coarsening on long time 
scales as conditions are 'fluxional' [|29ll. The final 'dried-in' structures depend on the pathway of 
evaporative dewetting (i.e. on kT, /i, the interaction constants e^-, nanoparticle concentration cj) 
and mobility M) and range from labyrinthine to polygonal network structures or holes in a dense 
particle layer. Typical snapshots from the evolution of a layer of solution are shown in Fig. |7J As 
before, the evaporative dewetting process of the solvent follows either a nucleation and growth 
[Fig.[7Jb)] or a spinodal-like [Fig.jT^a)] pathway. Following the mean field argument as above we 



14 




FIG. 7: Snapshots from typical evaporative dewetting processes with nanoparticles in the case of (a) a 
spinodal-like process at /i = —2.55 and (b) of nucleation and growth of holes at /i = —2.3. Starting 
condition homogeneous liquid film with homogeneously distributed particles. The number of Monte Carlo 
steps is indicated in the individual panels (top: 10, 15, 20, 100 bottom: 100, 200, 300, 400); e nn = 2, 
e n i = 1.5, M = 50, <j) = 0.2. The remaining parameters are as in Fig.[6j Particles are black, liquid is gray 
(green online) and the empty substrate is white. 

expect for the limiting curve fi ms (kT) separating the two processes in the parameter plane again 
a shift by — 2{e n i — l)cj). However, due to fluctuations this border is elusive when scrutinized 
numerically, therefore we will not try to compare simulations to the prediction. 

At first sight, however, one might get the impression that the particles act as a type of passive 
tracer that preserves the transient volatile dewetting structures of the solvent. This idea was put 
forward in Refs. [[37l |38l [391 in the context of experiments on dewetting aqueous solutions of 
macromolecules and provides an explanation for some of the basic features of the observed net- 
work structures. One can also employ this hypothesis to explain some of the structures observed 
in the present study, such as the network and spinodal structures shown in Fig. [7} The simulations 
indicate, however, that the nanoparticles are not simply passive tracers. Although the particles 
primarily just follow the solvent, they play an important role in several phases of the process. 

First, the particles may influence the nucleation process in the metastable parameter region, 
beyond the 'shift' in fi predicted by the mean field considerations above. Although our simulations 
seem to indicate that the nucleation rates actually depend on the particles themselves, we cannot 
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FIG. 8: Snapshots illustrating the unstable growth of a nucleated hole in a film of nanoparticle solution. The 
strong branched fingering caused by the nanoparticles can be seen clearly. From top left to bottom right the 
number of Monte Carlo steps is 450, 1200, 1800, and 3000. Nanoparticle related parameters are e nn = 2, 
e n i = 1.5, = —2.3, M = 5, (p = 0.2. The remaining parameters are as in Fig. |6j Particles are black, 
liquid is gray (green online) and the empty substrate is white. 

quantify the effect at present. We leave this question for future investigations and focus here on 
the second and more pronounced influence. Second, for a low mobility of the particles, i.e., for a 
slow diffusion of particles as compared to the evaporation of the solvent, the nucleated holes might 
grow in an unstable manner as illustrated in Fig.[8| The transverse instability of the dewetting front 
responsible for this effect is analyzed in more detail in the following section. 

V. THE FINGERING INSTABILITY 

A straight evaporative dewetting front without nanoparticles present moves with a constant 
mean velocity v (averaged over the transverse direction). The velocity only depends on the driving 
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chemical potential. The local velocity, however, fluctuates with temperature. The fluctuations are 
illustrated in the column of snapshots at the right of Fig. [5j Therefore, the local front position 
also fluctuates around the steadily advancing or receding mean position. The velocity is easily 
measured numerically for driving forces small enough that the nucleation of additional holes is 
not very probable during the time of measurement. Here this is the case for /i > —2.4. The 



dependency of v on /i is given for comparison below in Fig. [T2J One finds that the velocity 
increases linearly with the driving force. 

Adding diffusing nanoparticles to the solvent changes the front behavior dramatically as it may 
become transversally unstable. Such an unstable front does not only fluctuate around a contin- 
uously receding mean position but fingers stay behind permanently (for the front instability of a 
growing circular hole see Fig. [8] above). To elucidate the underlying mechanism we will in the 
next paragraphs discuss the influence of the individual system parameters on the instability, i.e., 
the influence of chemical potential /x, particle concentration 0, and mobility M. 

For each dependency we present snapshots of the final dried in nanoparticle structure corre- 
sponding to the dried in patterns observed in experiment. We add, however, thin lines that cor- 
respond to positions of the dewetting front at equidistant times throughout the evolution. The 
final branched patterns we characterize in a rather simple but effective manner. We determine an 
averaged finger number by dividing the side length of the computational square domain by the 
average finger number. The latter is obtained by counting the fingers on each line orthogonal to 
the mean direction of front motion (here the y direction) and averaging over the y-range where 
fingers exist. The underlying hypothesis that the finger patterns are 'stationary' will be discussed 
and checked below. Attempts to determine a mean distance of fingers employing a 1-D Fourier 
transform unfortunately did not give meaningful results because of the strong side-branching (see 
snapshots below). Note, that the finger number corresponds to a wave number rescaled by 2tt/L 
where L is the side length of the computational domain (that is fixed throughout this section and 
corresponds to about 1.2/xm.). 



A. Dependence on chemical potential 

First, we focus on the influence of the driving force, i.e., the chemical potential. Fixing the 
particle concentration to cp = 0.1 and the mobility to M = 20, the final structures observed for 
various values of fi are shown in Fig. [9] whereas Fig. [T0| shows the dependence of the mean finger 
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FIG. 9: Final 'dried-in' branched fingering structure for evaporative dewetting of a nanoparticle solution 
for different driving forces, i.e., chemical potential (a) /i = —2.1, (b) /i = —2.2, (c) fi = —2.3, and (d) 
fi = —2.4. Thin lines correspond to positions of the dewetting front at equidistant times with At = 333 
(a), and At = 166 (b, c, d). The domain size is 1200 x 1200, M = 20, (j) = 0.1, kT = 0.2, e nn = 2, and 

S n i 1.5. 

number on the chemical potential. The driving force increases for decreasing /i and the number 
of fingers increases roughly linearly. Again, for /x < —2.4 nucleation becomes more probable and 
we observe a random polygonal network. The side-branches of the polygonal pattern clearly result 
from the transverse instability [see Fig.[9]^d)]. 

Further inspection of Fig. [9] results in the following observations: 
(i) Both the wavenumber and the growth rate of the front instability increase with decreasing fi. 
A faster growth translates into longer fingers in the 'dried in' structure, i.e., fingers that extend 
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FIG. 10: Dependence of mean finger number < / > of the final dried in structures (determined as described 
in the main text) on chemical potential ji. Every data point corresponds to the mean value of 7 simulations; 
error bars indicate the corresponding standard deviation. The solid line corresponds to a linear mean-square 
fit with coefficients given in the plot. The remaining parameters are as in Fig.[9j 

further in Fig. [9| Starting from the same initial front position, the fingers appear much earlier in 
panel (c) than in panels (b) or (a). 

(ii) The front contours taken at equidistant times indicate that the front is already strongly disturbed 
before the instability becomes manifest in a deposited finger, i.e., parts of the front slow down 
before actually stopping (and depositing material). This is especially well visible in Fig. ^ a). 

(iii) The dried-in fingering structure is not entirely frozen as the chosen parameters values are in 
a 'fluxional' regime E9l . This is especially well visible at the tips of the fingers in Fig.[9fa) that 
continue to retract very slowly, or in Fig.|9fb) where a thin part of the rightmost finger breaks and 
slowly retracts as well. Eventually, this will lead to a long-time coarsening of the finger structure. 

(iv) It is remarkable that even when fingers are left behind the local distance between consecutive 
front contours seems to be constant in each of the panels [beside the regions discussed in (iii)]. 
One expects, however, that the front velocity depends on the local particle concentration and the 
mobility of the particles. The observation of a constant velocity indicates that the front instability 
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is related to an auto-optimization of the front velocity. This is analogous to an effect described for 
dewetting polymer films Bill . There, a moving dewetting front 'expels' part of the liquid rim and 
leaves the liquid behind in the form of deposited droplets. The instability avoids a slowing down 
of the dewetting front. Here, some of the particles that the evaporative dewetting front collects are 
expelled and deposited in fingers when too many are collected by the front. 

Measurements of the front velocity for different \i and particle concentrations indeed show that 
the mean velocity is constant and depends linearly on /i (see Fig. [T2| below). 



B. Dependence on particle concentration 

We have also examined the influence of the nanoparticle concentration (j) on the fingering. 
Fixing the chemical potential at ji = — 2.2 and the mobility at M = 10, the final structures are 



shown for various values of cp in Fig. [TTJ From an initial inspection of Fig. [TT] one might get the 
erroneous impression that the finger number increases with concentration. However, the statistics 
of many runs shows that the number is nearly independent of the particle concentration. The 
fingers only become thicker with increasing concentration. This indicates that the fingering is 
controlled by the 'dynamical' parameters responsible for the ratio of the time scales for diffusion 
of the particles and motion of the evaporative front (i.e., mobility and chemical potential). As the 
influence of the chemical potential has already been established in Section [VA] we focus below in 



Section V C on the influence of mobility. 

Before we do so we summarize in Fig.[l2]results for the dependence of the global average front 
velocity on chemical potential and particle concentration. For all investigated concentrations the 
velocity decreases linearly with increasing chemical potential, i.e., with decreasing driving force. 
If we fix the chemical potential, the dewetting process is slower for higher particle concentrations. 
This implies that the finger number is not directly related to the mean front velocity, a result that 
will be discussed in the Conclusion. 



C. Dependence on mobility 

Fixing the particle concentration at cj) = 0.1 and the chemical potential at fi = —2.2, we show 



in Fig. U3\ the final structures observed for selected values of mobility, M. Fig. 14 shows the 



dependence of mean finger number on mobility. One sees at once that mobility is very influential 



20 



FIG. 1 1 : Final 'dried in' branched fingering structure for evaporative dewetting of a nanoparticle solution 
for different nanoparticle concentration (a) (p = 0.05, (b) <p — 0.10, (c) (p = 0.20, and (d) cp = 0.30. Thin 
lines correspond to positions of the dewetting front at equidistant times with At = 166 (a, b), At = 333 
(c) and At = 500 (d). The domain size is 1200 x 1200, M = 10, /i = -2.2, kT = 0.2, e nn = 2, and 

S n i 1.5. 

in the fingering process. For decreasing mobility of the particles the number of fingers increases 
strongly. As was observed for the influence of the chemical potential, the growth rate of the 
instability increases with increasing finger number, i.e., decreasing mobility. This is manifest in 
Fig.[l3]by the later emergence of fingers for larger M, i.e., in the increasing distance of the finger 
tips from the upper domain boundary. 

The systematic picture that emerges from the preceding three sections allows us to identify the 
primary factors influencing the fingering instability. It becomes obvious that they are of dynam- 
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FIG. 12: Dependence of mean front velocity on chemical potential fi and for selected particle concentrations 
(j) as given in the legend. The thin lines correspond to linear regression fits with coefficients as given in the 
plot. The used domain size is 1200 x 1200, M = 20, kT = 0.2, e nn = 2, and e n \ = 1.5. The velocity is 
averaged over 7000 MC steps or the number of steps it takes to 'dry ' the substrate, whichever is smaller. 

ical nature. Mobility is the major player. The faster the nanoparticles can diffuse away from the 
dewetting front, the less likely the instability becomes as the front is not fast enough to 'collect' 
the particles efficiently. This directly corresponds to the observation that the fingering becomes 
more intense with increasing driving force, i.e., lower chemical potential. Then the relative dif- 
fusivity is hold constant (mobility M) but the mean front speed increases with decreasing fi. The 
greater the number of particles collected at the moving front, the more likely become lateral den- 
sity fluctuations that lead to self-amplifying variations of the local front velocity. Note, however, 
that up to this point we kept all the interaction parameters fixed. We expect them to influence 
the feedback loop because, for example, a stronger particle-particle interaction e nn should support 
clustering. But one has to be cautious because the process is of dynamic character. Therefore 
it is difficult to predict whether a larger e nn implies less or more fingers. This question will be 
discussed below in Section IVEl 

Next, however, we will check the assumption of stationarity of the finger patterns that we have 
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FIG. 13: Final dried in branched fingering structure for evaporative dewetting of a nanoparticle solution 
for different nanoparticle mobilities (a) M = 2, (b) M = 5, (c) M = 10, and (d) M = 20. Thin lines 
correspond to positions of the dewetting front at equidistant times with At = 166. The domain size is 
1200 x 1200, (f) = 0.1, ii = -2.2, kT = 0.2, e nn = 2, and e nl = 1.5. 

used above to introduce the measure of the mean finger number. 



D. Stationarity of finger pattern 



In Sections |V A| to |VC| we have used the mean finger number to characterize the parameter 
dependency of the fingering patterns. We noted that it only represents a trustworthy measure if the 
fingering process is stationary. By "stationary" we mean that the main properties of the pattern 
do not depend on the streamwise position. In other words the moving dewetting front moves on 
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FIG. 14: Dependence of mean finger number < / > of the final dried in structures on particle mobility 
M. Every data point corresponds to the mean value of 7 simulations; error bars indicate the corresponding 



standard deviation. The remaining parameters are as in Fig. 13 



average with a constant velocity and deposits on average the same amount of material in the same 
number of fingers. 



To check this strong assumption we plot in Figs. [15] and [16] the dependence of finger number 
on the streamwise co-ordinate (y-co-ordinate) for various values of the mobility and chemical 
potential, respectively. To smooth out fluctuations each curve represents the mean of 28 runs. As 
the fingers stay fixed after being deposited behind the front, an advance in y also indicates passing 
time. We note in passing that all 'snapshots' above that show dried-in final structures can as well be 
read as space-time plots tracking the particle-deposition at the moving dewetting front. Inspecting 
the curves we note that they have some properties in common: (i) finger number starts from zero 
at some value y = y s well behind the initial position of the front at y = 60; (ii) the finger number 
grows exponentially on a typical length scale related to the (linear) growth rate of the fingering 
instability; (iii) a maximal finger number is reached before the number decreases again and settles 
onto a stationary level. The decay corresponds to a small coarsening of fingers and can be seen 
in several of the snapshots above, (iv) Fluctuations around the stationary level are small, i.e., on 
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FIG. 15: Dependence of finger number on streamwise coordinate (i.e., y-direction) for various values of 
mobility as given in the legend. As the structures are fixed when they stay behind the front increasing y 
also indicates increasing time. The figure may be read as a space-time plot of the fingering at the moving 
dewetting front. Each curve corresponds to the average of 28 runs. The remaining parameters are as in 



Fig. 13 



average as many new fingers are created (branch tips in the final image) as vanish when fingers 
join pairwise (at the branching points). 

Our conclusion is that the mean finger number as used above is a valid measure because the 
y-range where the initial overshooting occurs is relatively small, and also the overshoot itself is 
normally below 20% of the mean stationary value. The resulting error is small as compared to 
the natural variance that we found between different runs. Note, however, that for several curves 



in Fig. [16] the decay after the initial maximum is relatively slow and the stationary value is not 
reached for our system size. This corresponds to a very slow coarsening of fingers as the front 
recedes. 
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FIG. 16: Dependence of finger number on stream wise coordinate (i.e., ^/-direction) for various values of the 
chemical potential as given in the legend. Each curve corresponds to the average of 28 runs. The remaining 
parameters are as in Fig. [9] 

E. Influence of interaction parameters 



We have mentioned in Section [TV] the hypothesis (put forward in ll39l for a related system) 
that the ongoing pattern formation and the properties of the receding dewetting front mainly result 
from the dynamics of the liquid-gas phase transition. The nanoparticles were seen as passive 
tracers. This describes well the situation in the spinodal-regime (taking into account that the 
mean field parameters are slightly changed by the particles). We have then shown that this picture 
has to be amended because the particles have an important influence during crucial phases of 
dewetting. Most notably, they are solely responsible for the transverse instability of the dewetting 
front leading to the branched finger patterns analyzed in detail in the present Section [V] 

The analysis up to now has been performed by fixing the interaction parameters at con- 
venient values: the liquid-liquid interaction serves as an energy scale, i.e., su = 1; the particles 
attract each other slightly more than particles attract liquid (e nn = 2, ^ = 1.5). This ensures that 
particle-liquid phase separation does not generally occur in the bulk solution, but also that particle 



26 



clusters are likely to form at higher particle concentrations. This resembles the situation in the 
experiments and was also used in most computations in Refs. Il26ll29ll . 
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FIG. 17: Dependence of the mean finger number on the particle-particle interaction strength e nn . The 
regions marked (i) to (iv) are discussed in the main text. The remaining parameters are: domain size 
1200 x 1200, kT = 0.2, M = 20, fi = -2.2, (p = 0.1, e ni = 1.5. The insets give typical snapshots 
obtained in the four different regions. Particles are black, liquid is gray (green online) and the empty 
substrate is white. 

In the following we will Analise the influence of e nn and e n \ on the fingering instability. In 
particular, we are interested in a possible coupling of the front instability and the demixing of 
particles and liquid. In our scaling (en = 1), demixing will occur for large e nn (at fixed e n i) or 
small e n i (at fixed £ nn ), i.e., for large ratios e nn I e n \. 



Fig. [T7] presents results for the average number of fingers when changing e nn in the range 
between 1.4 and 2.5. Inspecting the figure one can distinguish four regions: 

(i) At small e nn < 1.95 the particular strength of the interaction has nearly no influence; only 
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FIG. 18: Dependence of the mean finger number on the particle-liquid interaction strength e n \. The regions 



marked (i) to (iv) are discussed in the main text. The remaining parameters are as in Fig. 17 and e nn = 2.0. 



a slight decrease of the finger number with increasing e nn is discernible. We call this the 
'transport regime' as only the transport properties play a role. 

(ii) In an intermediate region 1.95 < e nn < 2.15 the finger number increases steeply by about 
a factor of three. Fingers emerge at the receding front and there is no demixing at all going 
on in the liquid bulk behind the front. What we see is still a pure front instability that is, 
however, strongly influenced by the ratio of interaction parameters. We interpret this as 
a front-induced demixing that leads to fingering. In other words, the front itself acts as a 
'nucleation site' for demixing. We call this a 'demixing-induced front instability'. Note, 
however, that fingering still leads to a strongly anisotropic fingering pattern. Although the 
fingers break at some places they do in general not break up into small nanoparticle islands. 

(iii) Increasing e nn further above 2.15, the finger number decreases again till about < / >= 15 
at about e nn = 2.25. This is mainly a geometric effect resulting from our one-dimensional 
finger counting routine. It reflects the fact that the fingers break up increasingly and the 
'dried in' nanoparticle pattern starts to look more and more isotropic. Demixing of particles 
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and liquid occurs already in the the bulk liquid behind the front. 

(iv) Beyond e nn = 2.15 fluid and particles demix already in the bulk in a homogeneous man- 
ner, the remaining dried-in structure of nanoparticle islands is isotropic, and the emerging 
demixing length scale decreases with increasing e nn . The 'finger number' is one way to 
count the islands. Note, however, that it is not an adequate measure to study those isotropic 
structures (that are not of central interest here). The front resembles a liquid front receding 
inside a porous medium (formed by the nanoparticle islands). 



Our interpretation of the ongoing physical processes is confirmed by Fig. [18] showing the de- 
pendence of the average number of fingers on the liquid-particle interaction strength e n \ in the 
range between 1.0 and 2.0 for fixed e nn = 2.0. A similar sequence of regions (i) to (iv) is found; 
this time however, with decreasing interaction parameter. This agrees well with the considerations 
at the start of this section. 



VI. CONCLUSIONS 

The present work has focused on pattern formation observed in various experimental settings 
involving dewetting and drying nanofluids on solid substrates. In addition to polygonal networks 
and spinodal-like structures, branched structures have been reported in the experiments. We have 
employed a kinetic Monte Carlo model to study pattern formation driven by the interplay of evap- 
orating solvent and diffusing nanoparticles. A justification of the usage of a model that only in- 
cludes dewetting by evaporation but not by convective motion of the solvent has been given based 
on experimental observations and scaling considerations derived from a mesoscopic continuum 
model. 

The model has first been used to analyze the influence of the nanoparticles on the basic dewet- 
ting behavior, i.e., on spinodal dewetting and also on dewetting by nucleation and growth of holes. 
It has been found that the 'classical' hypothesis that the solute mainly 'decorates' and 'conserves' 
the volatile dewetting structures of the solvent has to be amended in some important regions of 
parameter space. While it is true that the nanoparticles help to image the basic dewetting patterns 
(such as the labyrinthine structure resulting from spinodal evaporative dewetting or the random 
polygonal network resulting from nucleation and growth), they shift the bimodal and spinodal line 
in the chemical potential by a small amount that can be estimated using a mean field argument. In 
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consequence, they influence the nucleation rate. Although in the present work we have not focused 
on this aspect, it should be further scrutinized in the future. More importantly, the nanoparticles 
strongly destabilize straight or circular dewetting fronts in the nucleation regime. The main body 
of the paper has entirely focused on a numerical study of the underlying transverse front instability. 

We have analyzed the dependence of the characteristics of the branching patterns on the driving 
chemical potential, the mobility or diffusivity of the nanoparticles, and their concentration. As a 
result we have found that the mean number of fingers is almost independent of the particle con- 
centration. At a higher concentration one finds a very similar number of slightly thicker fingers. 
The most influential parameters are the particle diffusivity and the chemical potential. A decreas- 
ing chemical potential [mobility] leads to denser [less dense] finger patterns that develop earlier 
[later]. Note that a decreasing chemical potential corresponds to an increasing driving force, i.e., 
to a larger mean velocity of the dewetting front. As a result, we have drawn the conclusion that the 
crucial factor determining the instability characteristics is the ratio of the timescales of the differ- 
ent transport processes. In particular, a small ratio of the velocity of the dewetting front and the 
mean diffusion velocity of the particles renders the front more unstable. If the particle diffusivity 
is too low the particles are 'collected' at the front. In consequence, the front is slowed down, and 
an unstable stratification of concentration evolves that triggers a self-optimization of the front mo- 
tion: The front expels particles into fingers that stay behind. In this way it can maintain its velocity 
at a nearly constant value. A similar effect was described for dewetting polymer films Bill . There, 
the dewetting front expels liquid from the growing rim which collects the dewetted polymer and 
the liquid is left behind in the form of deposited droplets. The present system, however, allows 
one to also study the branching of the finger structure. 

Beside the dependencies on particle density, mobility and driving force we have also investi- 
gated the influence of the energetics of the system, i.e., the dependence of the fingering on the 
interaction strengths. We were especially interested in a possible coupling between front mo- 
tion and instability on the one hand and a liquid-particle demixing on the other hand. On phys- 
ical grounds we have distinguished three different regimes: at small particle-particle interaction 
strength or large particle-liquid interaction strength the fingering is nearly independent of the in- 
teraction constants. That is, it is a purely dynamic instability. We have called this the 'transport 
regime'. At large particle-particle interaction strength or small particle-liquid interaction strength 
demixing of liquid and particles occurs already in the bulk liquid. The retraction of the liquid front 
then resembles the dewetting of a two-dimensional porous medium. At intermediate interaction 
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strengths, the finger density strongly depends on the energetics of the system. We have found that 
particles and liquid demix at the moving front rendering it transversally unstable. The resulting 
length scale is determined by the dynamics and the energetics of the system. We have called this 
a 'demixing-induced front instability'. 

In the light of the numerical findings presented here the experimental results described above 
may be understood in terms of changes in nanoparticle mobility resulting from differences in the 
length of the carbon chains of the thiol ligands and the interpenetration of the ligands of different 
particles. Particles with short thiol chains diffuse faster and interpenetrate less. Therefore nearly 
no fingering is observed for C 5 and C 8 . Longer chains lead to slower diffusion implying better 
developed fingers for Ci and Ci 2 . It is likely that the decrease in fingering for Ci 4 -passivated 
nanoparticles results on the one hand from end gauche defects in Cu chains that tend to produce 
less pronounced interdigitation than for Ci 2 thiols [31 J. On the other hand longer chain length will 
also lead to greater core-core separation. This might drive e nn down and shift the system from 
a front demixing-induced front instability towards the transport regime where the finger number 
is lower. The present experiments, however, do not allow a direct comparison of finger numbers 
with the simulations as the radial geometry of the 'normal' nucleated holes does not facilitate a 
quantitative analysis. Future experiments that study fingering using imposed straight dewetting 
fronts would be of high interest. 

We had based the quantitative analysis of the fingering on the assumption of stationarity of the 
fingering process, i.e., on the assumption that the main properties of the pattern do not depend on 
the streamwise position. A detailed test of the assumption has shown that it is valid as the fingering 
pattern is stationary after an initial exponential growth in finger number that peaks at a maximal 
value before settling on a slightly lower level via finger coarsening. Following the transient phase 
the finger density fluctuates around a stationary state. This observation not only validates the 
quantification used but also confirms the hypothesis of auto-optimization of the front velocity: the 
front collects particles at a constant rate and leaves them behind in deposited fingers at the same 
rate. However, not only is the overall rate stationary but so too are the spatial distribution of fingers 
and its dynamics including the creation of new fingers (finger tips) and the annihilation of fingers 
by merging (branching points). 

We have checked that the results obtained are generic in the sense that small changes in the 
model set-up do not change the behavior qualitatively. In particular, we used an independently 
written code for 1 x 1 nanoparticles to ensure that the nanoparticle size is not a crucial element. The 
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present work assumed 3x3 nanoparticles to allow for simple comparison with already published 
results [[26l[27l[29l. To reach the same density of fingers for 1 x 1 nanoparticles the mobility has 
to be smaller than in the 3 x 3 case. 

A crucial rule used in the algorithm is the contrast in the mobility of the nanoparticles, i.e., 
the assumption that they do not diffuse onto the dry substrate. If one lifts this rule and tries 
to bound the nanoparticles stronger into the liquid by increasing the particle-liquid interaction 
strength, i.e., one strongly increases the energetic bias of the particles towards the retracting liquid, 
no front instability has been found. This explains why Ising-type models for grain growth that 
include mobile diffusing impurities, to our knowledge have never reported instabilities of domain 
boundaries caused by the diffusing impurities [[67l [68l |69]]. In all these models the impurities 
diffuse equally well in the two phases. They might, however, be energetically biased towards 
one of the phases or towards the phase boundary. Our present results on evaporating nanoparticle 
solutions suggest that unstable domain boundaries might also be found in the ordering dynamics 
of a ferromagnetic system with diffusing impurities under the influence of an external field if the 
mobility of the impurities strongly depends on the phase. 

Finally, we discus the limitations of the presently used kinetic Monte Carlo model that is based 
on the assumptions that the dominant processes can be captured in a two-dimensional setting ne- 
glecting the film thickness of the evaporating film and that the only relevant dynamics corresponds 
to particle diffusion and solvent evaporation. It was argued above, based on the general structure 
of a thin film model, that for film thicknesses < 10 nm, the convective motion of the solution 
can be neglected as compared to that due to evaporation. We have also laid out in which way re- 
cent experiments support this approximation. However, the experiments clearly show effects that 
cannot be described employing the present model. These effects are associated with the macro- 
scopic dewetting front which recedes before the ultrathin 'postcursor' film evaporatively dewets 
(and thus produces the different types of structure described here and elsewhere using variants of 
the two-dimensional kinetic Monte Carlo model [[26l|27l|29l[3T]|). The macroscopic front can also 
be unstable. That instability, however, obviously involves the convective motion of the solution. 
On the other hand one notices in the last three panels of Fig. [T] a bimodal network structure. The 
small scale network one is able to describe with the present model. The evolution of the large scale 
(mesoscopic) network, however, occurs at larger film thicknesses and involves convective motion 
of the solvent. Different models are necessary to describe those processes. Several mesoscale 
dynamic models such as , e.g., a dynamic density functional theory for nanoparticle and liquid 
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densities and thin film models are possible candidates for such models and will be investigated in 
the future. 
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